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Abstract 

A method is suggested allowing for the improvement of accuracy of self-similar factor 
and root approximants, constructed from asymptotic series. The method is based on 
performing a power transform of the given asymptotic series, with the power of this 
transformation being a control function. The latter is defined by a fixed-point condition, 
which improves the convergence of the sequence of the resulting approximants. The 
method makes it possible to extrapolate the behaviour of a function, given as an expansion 
over a small variable, to the region of the large values of this variable. Several examples 
illustrate the effectiveness of the method. 
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1 Introduction 



In the majority of realistic computational problems, the sought function, satisfying a very 
complicated set of equations, cannot be defined for the whole range of its variable, but can 
be found only for asymptotically small values of this variable. At the same time, the most 
interesting could be the behaviour of the function at very large values of the variable. This is 
the standard situation in extrapolation problems, repeatedly appearing in various applications. 

Suppose we are looking for a real function f(x) of a real variable x G [0, oo). By means 
of perturbation theory or an iterative procedure, we can find the behaviour of this function 
at asymptotically small x — > 0. But what of the most practical interest in many cases is the 
behaviour of f(x) at very large x, say as x — > oo. It is this the most difficult extrapolation 
problem that we address in the present paper: How, knowing the behaviour of f(x) only at 
x — > 0, to define the value /(oo) for x — > oo. 

There exist several extrapolation methods, among which the most known are the Pade sum- 
mation [1], Borel and Pade-Borel summations [2], and the optimized perturbation theory [3]. 
The latter was, first, advanced in Ref. [3] and nowadays is widely employed for various applica- 
tions, as can be inferred from the review works [4,5]. Another extrapolation method is based on 
the self-similar approximation theory [6-11]. Using the techniques of this theory, supplemented 
by the fractal transforms [12,13], we have recently derived a novel types of approximants al- 
lowing for an effective summation of power series and, hence, for their extrapolation. These 
are the self-similar exponential approximants [14,15], self-similar root approximants [16-18], 
and self-similar factor approximants [19,20]. The exponential approximants are appropriate for 
extrapolating functions corresponding to exponentially varying processes, while the root and 
factor approximants suit well for extrapolating functions with power-law behaviour. The self- 
similar approximants were shown to be simpler and more accurate than the Pade approximants 
[16-20]. 

In the present paper, we aim at improving further the accuracy of the factor and root approx- 
imants by introducing a control function through a power transform of the initial asymptotic 
series. The new technique is illustrated by several examples of extrapolation from f(x) at x — > 
to /(oo). 



2 Self-similar summation of power transforms 

Assume that the behaviour of the sought function f(x) is known only for asymptotically small 
x — > 0, when the function can be represented as an expansion in powers of x, being approximated 
by the series 

k 

fk{x) = Y, a nX H , (!) 
n=0 

where k — 0,1,2, . . .. Without the loss of generality, we may consider such expansions for which 
a = 1, so that /(0) = 1. Really, in the case when 

f(x) - fo(x) (l + aix + a 2 x 2 + . . .) , 

with a nontrivial fo(x), not expandable in a power series, we may always define 

7(x) = ^ 

/( j -/o(x)' 
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after which the series for f(x) acquire the form (1), where ao = 1. Our aim is, being based on 
the behaviour of f(x) at x — > 0, where it is approximated by the series (1), to find /(oo) at 
x — > oo. 

The novel trick, we advocate in this paper, is to introduce a control function m = m k {x) by 
means of the power transform 

P k (x,m) = fr(x). (2) 
Taking the power m of series (1), we reexpand the result in x obtaining 

P k {x,m) = j^b n {m)x n , (3) 

n=0 

with b n (m) defined through a n . A particular case of transform (2) is an inverted series with 
m — — 1, which we have considered earlier. However, fixing the power m is not the best choice 
and here we shall advance a more general and rigorous way of selecting m. Expansion (3) 
serves as a basis for constructing in the standard ways [12-20] the self-similar factor and root 
approximants. Even-order factor approximants [19,20] are defined as 

F 2k (x,m) = '[[(1 + A i x) n < (4) 
i=i 

and odd factor approximants can be represented as 

k 

F 2k+ i(x,m) = l + MlI(l + ^r , (5) 

i=i 

with the parameters A { = Ai{m) and rij = rij(m) defined by the accuracy-through-order proce- 
dure with respect to series (3). This means that Eqs. (4) or (5) are to be expanded in powers 
of x, and these expansions have to be compared with Eq. (3), equating the terms of like orders. 
For the root approximants [16-18], we have 

R 2k (x, m) = ((... (1 + Aix) ni + A 2 x 2 Y + ... + A k x k ) nk (6) 
in even orders, and 

R2k+i(x, m) = 1 + hx ((. . . (1 + A lX ) ni + A 2 x 2 ) H2 + ... + A k x k ) nk (7) 

in odd orders. The parameters Ai = Ai{m) and n» = rij(m) could be defined in two ways. If the 
behaviour f(x) at x — > oo would be known, this could be used for uniquely defining all parame- 
ters [4] . Another way is to determine these parameters by means of the accuracy-through-order 
procedure. The second way may yield multiple solutions, which, however, are usually close to 
each other [21]. In what follows, we shall present only the most accurate approximant. The 
advantage of the root approximants is their ability to catch rather complicated asymptotic 
behaviour at large x, including corrections to the main scaling. 

Defining the parameters Ai(m) and rij(m) by the accuracy-through-order procedure, we 
assume that the limit /(oo) exists and finite, which imposes an additional constraint on the 
sum J2 k =i rii- 
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After constructing a factor approximant Fk(x,m) or a root approximant Rk{x,m) for the 
power transform (3), we have to accomplish the transformation inverse to Eq. (2), thus, ob- 
taining either 

f k (x,m) = [F k (x,m)] 1/m (8) 

or 

r k (x,m) = [R k (x,m)} 1/m . (9) 

The improvement of the accuracy, as compared to the factor and root approximants not 
involving the power transformation (2), is achieved by defining a control function m = m k (x) 
from a fixed-point condition. In general, there exist several types of such fixed-point conditions, 
which, actually, are equivalent to each other [4]. Here we use the simplest of them, the minimal 
sensitivity condition, which gives either 

d 

— f k (x, m) = , m = m k (x) , (10) 

or 

d 

— r k (x, m) = , m = m k (x) , (11) 

depending on whether the factor approximant (8) or root approximant (9) is considered. With 
the so defined control function m k (x), we get either 

r k (x) = f k (x : m k (x)) (12) 



or 



r* k (x) = r k (x,m k (x)) . (13) 



From here, keeping in mind our main aim to extrapolate the sought function to the limit 
x — > oo, we obtain either 

^(oc)= hm f* k (x) (14) 

or 

r *k(°°) = J™, r K x ) ■ ( 15 ) 

Below, we shall illustrate the method by several examples, confronting the found approximants 
f k (oo) and r£(oo) with known values of /(oo). Note that for defining the limits (14) or (15), 
we, actually, do not need to have the whole function m k (x), but what we need to have is just 
a limiting value m k = m k (oo), which is a constant. 

3 Stirling series for factorial function 

Let us consider the factorial function 

f(x) = 4= JI'xV'T (l + -) , 

V27T V xj 

where T(-) is a gamma function. As x tends to zero, one has 

fix) ~ (x - 0) . 
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Therefore, we define the reduced function 

fix) = ^fx f(x) , 

whose small-a; expansion has the form of Eq. (1), so that, as x — > 0, then 

fix) ~ 1 + a±x + a 2 x 2 + a 3 x 3 + a 4 x 4 + a 5 x 5 . 

The expansion coefficients are 

1 1 139 571 163879 

ai - — , a 2 - — , a 3 - - 51g40 , a 4 - - 2488320 > a 5 - 2 09018880 ' 

Thence, we shall apply the procedure of Section 2 to the function fix), and at the end, we will 
return to the sought function f(x) = f(x)/y/x, looking for the limit /(oo). The exact limit for 
the factorial is 

/(oo) = -= = 0.398942. 

V 27T 

Following the method, described in Sec. 2, we find 7-5(00) = 0.458, whose error, as compared 
with the exact /(oo), is 15%. The factor approximant f£{oo) = 0.406 is much better, with 
an error of only 2%. Comparing this with the Pade approximants, we should remember that 
these are not uniquely defined, yielding for each given order a whole table of approximants [1]. 
One often considers solely the diagonal approximants. For the present example, the diagonal 
Pade approximant -Pp/2] describes [22] the factorial-function limit /(oo) with an error of 14%. 
Thus, the factor approximant /^(oo) = 0.406 is the most accurate. It is worth emphasizing 
that a direct application of the factor-approximation technique, without involving the power 
transformation (2), would give the limiting value 0.169, which is a very bad approximation. 
Hence, employing the power transformation is a crucial point in improving the accuracy of the 
approximants. 



4 Debye-Hiickel function for strong electrolytes 

The function 

= \ ~ J* (1 - e~*) 

arises in the Debye-Hiickel theory of strong electrolytes [23]. At small x — > 0, this function 
possesses an expansion of the type (1), with 

11 11 1 



a l — ~ o a 2 — — , 03 — — — ; a A — 7777: , — 



3 12 60 360 2520 

We shall be interested in finding the limiting value /(oo) of the reduced function 

fix) = xDix) , 

whose exact limit is /(oo) = 2. 

Using the technique of Section 2, we get for the best root approximant t^oo) = 1.993, 
whose error is —0.4%. For the uniquely defined factor approximant, we find /^(oo) = 1.779, 
with an error of —11%. Note that without invoking the power transform, there are no real 
solutions for the sought limit. Thus, the usage of the power transformation (2) is principal 
here. The best Pade approximant, employing the same coefficients a n , gives the limit /(oo) 
with an error of —33%. 
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5 Critical temperature of Bose gas 



Bose-Einstein condensation in dilute Bose gas has attracted much attention in recent years 
(see reviews [24-27]). One of the interesting problems, which has been intensively studied, is 
the influence of atomic interactions on the shift of the critical temperature. One considers the 
relative variation of the critical temperature 

^JL = Il | 

due to weak atomic interactions, as compared to the condensation temperature 



2irh 2 

T n 



mki 



i 2 /3 



C(3/2) 



of the ideal homogeneous Bose gas. The lowest term in the expansion of the critical-temperature 
shift with respect to the small gas parameter 

7 = P 1/3 a s , 

where p is particle density and a s , scattering length, is commonly represented as 

AT C 

- ci7 (7^0). 



T 



c 



The coefficient c\ has been calculated by a number of various methods. Review of the related 
literature up to 2004 can be found in Refs. [25,27]. The most accurate are the results for c\ 
obtained by means of the Monte Carlo simulations and using the optimized perturbation theory. 
Less accurate are the results based on a renormalization-group approach [28,29]. Lattice Monte 
Carlo simulations by Arnold and Moore [30,31] give c\ = 1.32 ± 0.02 and by Kashurnikov et 
al. [32,33], Ci = 1.29 ± 0.05. Path integral Monte Carlo simulations by Nho and Landau [34] 
give Ci = 1.32 ± 0.14. A variant [5] of optimized perturbation theory, employed by Kastening 
[35-37], yields c\ = 1.27 ± 0.11, and the optimized perturbation theory used by Kneur et al. 
[38,39], results in c\ = 1.30 ± 0.03. Here, we shall calculate the coefficient c\ by means of the 
technique of Section 2. 

The coefficient C\ can be expressed as an asymptotic expansion 

ci(g) ^ ai9 + a 2 g 2 + a 3 g 3 + a 4 g 4 + a 5 g 5 

in powers of an effective coupling parameter [36], where 

ai = 0.223286 , a 2 = -0.0661032 , a 3 = 0.026446 , 

a 4 = -0.0129177 , a 5 = 0.00729073 . 
This expansion is valid for g — > 0. But the sought value of C\ is given by the limit 

ci = lim d(g) . 

g^oo 

Employing the factor approximants, complimented by the power transformation (2), we have 
/ 5 * (oo) = ci = 1.09, which is close to the values found by other methods. Summing the strong- 
coupling expansion with the help of the root approximants and defining the parameters from 
the weak-coupling expansion, we get rj^oo) = 1.19. 



6 



6 Structure factor of branched polymers 

The structure factor of three-dimensional branched polymers is given [40,41] by the confluent 
hypergeometric function 

s(«) = fi(4§«) . 

We shall consider the reduced function 

f(x) = xS(x) , 

whose limit /(oo) = 1/3 is finite. At asymptotically small x — > 0, the function f(x) possesses 
an expansion of the form (1), with a = 1. Several other expansion coefficients are 

ai = -1 , a 2 = 0.6 , a 3 = -0.257 , a 4 = 0.086 . 

The best approximant, obtained by the method of Section 2, is /|(oo) = 0.329, whose error is 
—1.3%. This is much more accurate than the best Pade approximant of the same order having 
an error of —266%. 



7 Discussion 

In this paper, we suggested a method for improving the accuracy of self-similar approximants 
by introducing a control function through the power transformation (2). As is shown by several 
examples, the accuracy really becomes essentially better. 

Here we have concentrated our attention on the extreme extrapolation problem, when from 
the behaviour of a function f(x) at asymptotically small x — > one has to find the limit 
/(oo) for x — > oo. This extrapolation problem is one of the most difficult. If we are able to 
accurately predict the behaviour of a function f(x) at x — > oo, then, as is clear, it is even easier 
to approximate its behaviour for finite x. 

As an illustration of the latter statement, we may consider the expansion factor of a polymer. 
The properties of polymers are of great importance for a variety of applications [42]. Let us, 
for example, consider the expansion factor a(z) for a three-dimensional polymer chain with 
excluded- volume interaction, where z is a dimensionless coupling parameter [43,44]. From an 
asymptotic series of the type (1), derived by means of perturbation theory [43], we construct 
the expansion factor 



a*(z) = 1.5286^ 



0.383 

0.3543 ' 



/ , \ -0.0749 n noKO 

(l + 0.15522- 1 ) • 0.3302: 1 



This is obtained by considering a large- ,2 expansion, resumming it by means of the root ap- 
proximants, and determining all unknown exponents and amplitudes from the weak-coupling 
expansion of fourth order. The strong-coupling exponent 

1 1 , In a(z) 
u=- + - hm - 

2 4 z^oo i n z 

for our approximant a*(z) is v — 0.5886, which coincides with the value found numerically 
[43-45]. That is, the approximant a*(z) possesses a correct scaling behaviour. It also gives a 
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nontrivial correction to the scaling, with an exponent of —0.9552. The expression a*(z) is valid 
for all z G [0, oo), differing from the known numerical values [44] by not more than 0.3%. 
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